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FOREWORD 

This report is one of a series of technical documents prepared for 
the Grand River Basin Water Management Study. The project described 
herein was undertaken by the Grand River Hydrology Sub-Committee at 
the request of the Grand River Implementation Committee. 

The material contained in these reports is primarily technical 
support information and in itself does not necessarily constitute 
policy or management practices. Interpretation and evaluation of 
the data and findings, in most cases, cannot be based solely on this 
report but should be analysed in light of other reports produced 
within the comprehensive framework of the overall study. Questions 
with respect to the contents of this report should be directed to 
the Co-ordinator of the Grand River Basin Water Management Study, 
Water Resources Branch, Ministry of the Environment, 135 St. Clair 
Avenue West, Toronto. 
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SUMMARY 

1. The three watershed models M0EHYDR2, SQGRID and TVA are reviewed and 
evaluated for their applicability in the Grand River basin. 
Qualitative evaluations are given in Chapter 3, whereas Chapters 4, 5 
and 6 deal with calibration, prediction and discussion of results. 
Description of parameters involved in each model and sample outputs 
are presented in Appendices A, B and C. 

2. The models are calibrated on test basins: Eramosa River, Lutteral 
Creek and Nith River. Only the Eramosa River sub-basin serves as a 
common basin among all three models. The M0EHYDR2 and SQGRID models 
are calibrated and compared using the Lutteral Creek sub-basin and 
the former is calibrated and evaluated further in the Nith River 
sub-basin. 

3. A direct comparison between the M0EHYDR2 and the SQGRID models for 
hourly simulations shows the models yield results which are quite 
comparable. The SQGRID model predicts peak flows for specific storms 
with reasonable accuracy. The M0EHYDR2 is more attractive because 
less CPU time is required. For the Eramosa River sub-basin, the cost 
of simulation by the SQGRID model is approximately seven times that 
of the M0EHYDR2 model. 

4. On a seasonal basis, the daily flows predicted by the M0EHYDR2 model 
are satisfactory (See figure 5.1). The SQGRID model cannot be used 
for direct, daily simulations because the time step is restricted by 
the grid size chosen. 

5. The M0EHYDR2 model is recommended for daily flow simulations for 
periods up to one year and for situations where a continuous 
simulation of both daily and hourly flows is required in 
combination. In the SQGRID model, the computations are performed 
hourly in each of the grids in sequence and the computer costs are 
high. This model, therefore, is not suitable when the simulations 
are required over a long period. 

6. By interfacing with a suitable precipitation forecasting system, 
SQGRID should form a good flood forecasting model. The SQGRID model 
is also useful in studying the effects of changes in the distribution 
of some physical characteristics of the watershed, such as 
infiltration rates or land use. 

7. The TVA model is recommended for long term simulations where average 
monthly flows will satisfy the flow requirements. 

8. The qualitative aspects of the three models and their performance in 
simulating flows are summarized in the table on the next page. 
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MODEL QUALITIES AND SIMULATIONS OF FLOWS 

MOEHYDR2 SQGRID 
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1. Water balance studies 
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in watershed characteristics 
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7. Input data availability 
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1 . INTRODUCTION 

1 . 1 Purpose of Report 

Based on the need for integrated water management actions, the Grand 
River Implementation Committee is co-ordinating a study to provide 
the necessary information for planning and development of the water 
and related land resources in the ,Grand River basin. The study 
deals with the important aspects such as flood control, water 
quality, pollution control, waste assimilation, water supply and 
flow augmentation within a general framework of a water resources 
management program for the basin. Implementation and use of 
hydrological models for the simulation of hydrologic regimes and 
applications in water resources management form a part of the study. 

During 1976/77, three watershed models were implemented for 
comparison and evaluation of their applicability in the Grand River 
water management program. This report presents the results of that 
study. Subsequent to this study, another model, the National 
Weather Service River Forecasting System (NWSRFS) was implemented, 
the results of which shall be presented in a separate report. 

1.2 Watershed Models 

A watershed model is a streamflow simulation model in which due 
consideration is given to the physical processes taking place in the 
rainfall-runoff phase of the hydrologic cycle. In this case, the 
inputs are usually the meteorological variables and the model 
transforms the inputs into streamflows. 

The analysis of rainfall-runoff processes of a drainage basin is 
complex and simulation models require the development of 
relationships necessary for the conversion of meteorological inputs 
into runoff. Watershed models are normally based on a combination 
of deterministic and parametric techniques. 



2. REVIEW OF THE SELECTED STREAMFLQW SIMULATION MODELS 

2. 1 Selection of the Models 

Different kinds of water resources and hydrologic models are being 
developed by various agencies. These models are developed to aid in 
the solution of either a specific problem or a set of problem 
situations. The choice of a particular model for application in 
solving a specific problem may depend on the overall development 
goal and objectives of the project. The prerequisite to the choice 
of a hydrologic model requires a clear definition of the main 
objective of the water resources management program. Furthermore, 
the model selected should be sufficiently flexible in operation so 
as to allow for its use for several purposes, with minor 
modifications. 

The choice of a suitable watershed model for application in 
hydrological studies in the Grand River Basin is important. Though 
numerous watershed models are reported in literature (WMO, 1975) not 
all of them meet the needs and requirements of the hydrological 
investigation for the Grand River Basin Water Management Study. The 
watershed model, in addition to the simulation of river flows, 
should be capable of simulating, with reasonable accuracy, various 
other components of the hydrologic cycle so that the model can be 
used in ground water and other water balance studies. 

In this preliminary report, three models are investigated for their 
applicability in the Grand River basin. The choice of these three 
models is based on easy accessibility and familiarity of the models 
to the working group, and the computer hardware requirements. 
Furthermore, two of the models selected (M0EHYDR2 1 and SQGRID 2 ) 
were developed in Canada and were cited to give satisfactory results 
for basins in Canada. The third model (TVAj) was developed in the 
United States of America and has been applied for water resources 
management studies of small basins in the Tennessee Valley area. It 
can be noted that each of the models under consideration belongs to 
a different category; namely they can be considered as a lumped 
model, distributed model and an index model, respectively. 



1. Comprehensive watershed model (lumped, deterministic) - Ontario 
Ministry of the Environment 

2. Rainfall-runoff watershed model (distributed, parametric) - 
University of Waterloo. 

3. Total rainfall-runoff watershed model (lumped, index parametric) 
- Tennessee Valley Authority, Tennessee, U.S.A. 



2.2 Review of the Models 

2.2.1 Lumped Model-M0EHYDR2 

The model "M0EHYDR2" is a computer program package for the 
simulation of a comprehensive watershed hydrologic system. The 
model was developed by the Water Modelling Section of the Ministry 
of the Environment. The model is formulated to represent, as well 
as practicable, concepts of the interrelationships of the various 
hydrologic processes in the watershed (Logan, 1972). Test 
applications were completed in International Hydrological Decade 
(IHD) representative basins in southern Ontario. 

(a) Model Structure 

The schematized form of the model structure is shown in Figure 2.1 
As indicated in the figure, the precipitation (including excess 
snowmelt) is processed through the various functions of 
interception, surface depression and detention storages, 
infiltration, ground water and surface runoff components. The 
infiltrated portion of the precipitation, after satisfying the 
maximum storage requirements of the aeration and interflow zones, 
contributes to the accretion of the ground water storage. The total 
of surface flow, interflow and the baseflow is routed to produce the 
total runoff at the outlet of the basin. The model package 
comprises three major sub-models: one for daily simulation (LTSM), 
another for hourly simulation (STSM) during specified storm event 
periods and the other for simulating snow accumulation and melt 
(SNMELT). All of the three submodels are interfaced for continuous 
daily simulation for a period of up to 366 days. In the hourly 
simulation mode, any number of days corresponding to a selected 
storm event may be included as desired. The M0EHYDR2 model is a 
lumped, deterministic, parametric model; and the simulations are 
carried out for average input conditions of the basin. Estimates of 
the most sensitive model parameters are calculated based on a 
trial-and-error basis supplemented by automatic optimization 
procedures. In the automatic procedure, the function which is 
optimized is the sum of squares of the differences between simulated 
and observed flows. The optimization is carried out by an automatic 
"pattern search" optimization routine. 

(b) Sub-model LTSM 

The LTSM sub-model is used in the daily flow simulation mode of 
M0EHYDR2. In this mode, seven parameters are optimized, five of 
which are related to infiltration and surface runoff components, and 
the other two deal with ground water storage and baseflow 
calculations. A complete list and description of model parameters 
and constants used in the model are given in Appendix A. 
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(c) Sub-model STSM 

The STSM sub-model is used in the hourly flow simulation mode of 
M0EHYDR2. The method of computation is essentially the same as that 
contained in the LTSM sub-model, except that some different 
algorithms are used for some of the calculations which are performed 
hourly during the designated storm dates. The optimization routine 
considers five parameters, three of which are related to 
infiltration and the other two are used in the surface runoff 
routing calculations (Appendix A). 

(d) Sub-model SNMELT 

This sub-model simulates the snow accumulations during the winter 
months and calcuates the excess melt from the accumulated snow on a 
daily basis. The method of melt calculations is based on energy 
balance principles and the sub-model is programmed to function on 
its own. It is then possible to interface this sub-model with any 
comprehensive watershed model. Further details on the snowmelt 
routine are presented by Logan (1976). 

(e) Input Requirements 

The input requirements for model initialization and application fall 
into five different categories of specification. 

(i) Model Control Variables 

Model control variables include the types of data sets used, 
start and end date of simulation period (accommodating leap and 
non-leap year), specification for selected input storm event 
sequences and monthly variation of average daily maximum 
sunshine hours. 

(ii) Data Input Variables 

Data input variables are meteorological and streamflow data for 
use in calibration or only meteorological input for use in 
simulation and prediction. 

For the LTSM sub-model the input are daily precipitation (total) 
and estimates of potential evaporation. 

For the SNMELT sub-model, in addition to daily precipitation 
(rain and/or snow), daily input of air temperature, relative 
humidity, net radiation and wind speed are required. 

(iii) Optimization Control Variables 

These variables include the number of parameters to be 
optimized, number of observations, number of iterations, type of 
objective function to be used, incrementor, decrementor and 
limit values for the parameters and test level of convergence. 
They specify control for optimization of LTSM or STSM sub-model 
parameters and specification for the desired output summary 
options. 



(iv) Watershed Model Constants 

These include physical basin parameters (area, flow path length, 
etc.), empirical constants and coefficients (soil moisture field 
capacity, wilting point, snowmelt coefficients, etc.). 

(v) Initialization Variables 

These include initial basin storage condition (soil, ground 
water, snowpack, flows, etc.) and initial values for model 
parameters. Definitions of the model constants required for the 
input are listed in Appendix A. 

(f) Data Processing 

The necessary meteorological data can be obtained in the form of 
punched cards or on tapes from the federal, Atmospheric Environment 
Service, Canada. Depending on the distribution of the 
meteorological stations within and around the basin, weighted 
average values may be used. As a complement to the application of 
M0EHYDR2, a computer program was developed to read the magnetic tape 
supplied by the AES, containing the meteorological data for the 
different stations. The abstracted data for a particular variable 
(precipitation or temperature) are used in computing weighted 
averaged values for the basin. These are provided as punched cards 
in the input format required by the M0EHYDR2 model. The different 
station weights are pre-specified inputs to the computer program. 

Estimates of values for the various constants and physically based 
parameters of the model can be obtained from field measurements and 
published information. Additional and supplementary statistical or 
empirical analyses are required to determine estimates of a number 
of less sensitive model coefficients. 

(g) Output 

In addition to the basic outputs of optimized parameter values and 
the predicted runoff, the output from the M0EHYDR2 model includes 
several options in the printout. An important feature of the output 
is the complete information on the water balance of the river basin 
printed out for every day of the period of simulation, including the 
snowmelt periods. For hourly simulations, the same information is 
provided on an hourly basis. There are ten different output options 
available. 

A list of the different printout options and samples of printouts 
and plots are given in Appendix A. 



2.2.2 Distributed Model - SQGRID 

The square grid model is a distributed parametric watershed model, 
developed at the University of Waterloo for the estimation of river 
flows from meteorological data (Gupta, 197*0. The model is 
developed as an extension to the distribution of the meteorological 
data on a grid system superimposed on a basin (Solomon et al 1968). 
Each grid cell is treated on an individual basis as a "lumped 
sub-element" within which the physical model parameters are assumed 
to remain constant. 

(a) Model Structure 

The hydrologic processes considered in this model, are very similar 
to those presented in Figure 2.1. However the calculations are 
performed separately for each grid cell and the flows, computed 
numerically on the basis of water balance and continuity equations, 
are routed through the elements in a pre-determined hierarchical 
order. 

The size of a grid in the system can vary from 1 km x 1 km to about 
10 km x 10 km, depending upon the size of the watershed. The time 
step for simulation is 1 hour. It is not possible to use time steps 
greater than one hour because instabilities may occur when the grid 
size is small and the velocity of flow within a grid cell is large. 
The optimizing function is the sum of squares of differences between 
observed and computed daily flows, the same as in the M0EHYDR2 
model. However, an automatic optimization routine "CLIMB" is used 
instead of a pattern search algorithm. 

In the grid model, the contribution to the river flows from ground 
water storage is not considered. The baseflow contributions 
resulting from individual storm events is assumed to be negligible, 
especially at upstream reaches of the stream. However, baseflows 
are significant when the simulation is carried out over several days 
over a number of grid elements integrated at the overall basin 
outlet. In its original form, the square grid model does not 
include a snowmelt routine and thus is not applicable during winter 
months. However, it may be possible to interface a suitable snowmelt 
subroutine with the main simulation model. 

The computer package program comprises basin files for the basin 
drainage configuration and flow network hierarchy, for land use and 
distribution, for soil types and distribution and for the 
meteorological input data. A schematic representation of the 
structure of the computer program and a description of the basin 
files are given in Appendix B. 

(b) Basin Files 

BASFI A - Establish hierarchical level of grid system in drainage 

basin. 
BASFI B - Overlay the digitized land use and soil data on the grid 

system on file. 
BASFI C - Delineation of a sub-basin above any point in the system. 
BASFI D - Generates digital and gray map, showing spatial 

distribution of data. 
BASFI E - Calculates distributed information for soil and land use. 



(c) Input Requirements 

The input requirements can be divided into four different categories 
as given below: 

(i) Model Control Variables 

These relate to specification for designated basin grid 
location, start and end of simulation period, and optimization 
criteria. 

(ii) Data Input Variables 

These include hourly precipitation and air temperature, hourly 
flows for use in calibration and basin data file inputs (see 
BASFI description). 

(iii) Optimization Control Variables 

These include the number of parameters to be optimized, number 
of iterations, initial value and upper bound for each parameter. 

(iv) Initialization Variables 

These are initial basin soil moisture and subsurface flow 
storages. 

(d) Data Processing 

The necessary meteorological data can be obtained from the 
Atmospheric Environment Service, Canada. A computer program 
"EXTRPREC" is used to extract the data and an interpolation program 
"INTERP" is employed to redistribute precipitation and temperature 
input to each grid element in accordance with station location and 
distance to grid element. Daily evaporation rates are calculated, 
as a function of air temperature, using an empirical equation 
(Turc's formula). 

The values for a number of model constants can be obtained from 
field measurements and published information. 

(e) Output 

The simulation output contains a table of hourly estimates of 
evaporation, soil moisture and ground water storage indices, 
variable overland and channel lag coefficients and the simulated 
surface and subsurface components of the runoff. Also included in 
this table are the hourly inputs of precipitation and temperature 
for the grid cell under consideration, as calculated by the "INTERP" 
subroutine. A summary of the simulated daily flows, plot of 
observed vs. simulated flows and the parameter values are given at 
the end of the optimization run. There is an additional output 
option to have the above information printed out at a maximum of 
four different locations in the basin. 

A sample plot of observed vs. simulated flows and a sample printout 
are included in Appendix B. 



2.2.3 Index Model - TVA 

The TVA watershed model was developed by the Tennessee Valley 
Authority, mainly to aid the planning and development of water 
resources in the tributary areas of the valley. The model was 
developed and tested in the Upper Bear Creek Experimental Project. 
This model is an extension of the NWS Coaxial (graphical) Antecedent 
Precipitation Index method to predict storm runoff from rainfall 
(Tennessee Valley Authority, 1972). 

(a) Model Structure 

The model comprises a series of conceptual hydrologic compartments 
interconnected according to an input-output cascading sequence. 

The computer package program consists of a number of simulation 
algorithms. Lumped retention index relationships for soil and 
ground water storages are used to determine the fractions of daily 
precipitation input to surface and subsurface flows. The routing of 
surface and subsurface flow increments through basin storages is 
based on empirical lag coefficients. The package program is 
completed with an automatic optimization routine and a plotting 
routine. Five primary model parameters undergo automatic 
optimization, and estimates of another five secondary parameters are 
made by conventional analytical methods. The optimization routine 
also provides statistics of goodness of fit for the model. The 
routines are programmed for continuous daily simulation for a period 
of up to three years. 

The original model does not include a snowmelt routine. Therefore, 
a simple snowmelt sub-model based on simple degree-day method was 
interfaced with the main program for simulation during winter months 
for its application in the Grand River basin. 

(b) Input Requirements 

The input requirements can be sub-divided into the following groups. 

(i) Data Input Variables 

The necessary meteorological input for the simulation are daily 
precipitation, air temperature and monthly estimates of 
potential evapotranspiration. Daily streamflows are required 
for the calibration of the model. 

(ii) Optimization Variables 

These include the number of parameters to be optimized, number 
of observations, number of iterations, initial values and 
incrementor and decrementor values and convergence limits for 
the parameters. 



(iii) Model Constants 

These are drainage area, seasonal coefficients for interception 
storage and growth index for crops, land use and distribution 
indices for soil moisture and ground water availability. 

(iv) Initialization Variables 

These include initial basin storage conditions (soil moisture, 
ground water, etc.). 

A list and description of the model parameters and constants are 
given in Appendix C. 

(c) Data Processing 

Daily precipitation (rainfall) values from a number of stations are 
processed according to the Theissen polygon weight to determine 
weighted average input for the basin. In addition, a "precipitation 
lagging correction" is made to account for the difference in times 
at which the rainfall and streamflow observations are made. 

The program has a "lead-in" period of three months so as to 
dissipate the effect of incorrect starting values for initial soil 
moisture storage. During the lead-in period the streamflow data are 
not required. 

(d) Model Output 

The output of the TVA model computer program lists, in a tabulated 

form, daily summaries of rainfall input, components of simulated 

flows, actual evapotranspiration and volume storage changes for 
surface moisture, soil moisture and groundwater. 

A plotting option is available for daily rainfall input, observed 
and simulated flows. A summary of the optimized model parameters 
and a set of goodness-of-fit statistics are also given. Samples of 
printout and plots are included in Appendix C. 
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3. MODEL EVALUATION CRITERIA 

An evaluation of different watershed simulation models for their 
applicability in water resources management studies in a particular 
river basin, should be based on both qualitative and quantitative 
criteria. A qualitative evaluation can be a comparison of the 
conceptual aspects of the model and their usefulness in a water 
management program. A quantitative evaluation can be done by way of 
model calibration and verification through simulation of flows. For 
example, the models may be calibrated on the same set of flow data 
and tested subsequently for their predictive capability with the use 
of another data set (WMO, 1975). Chapters 4 and 5 of this report 
deal with such an evaluation. 

3.1 Qualitative Evaluation 

In a qualitative evaluation of the concepts and functions of the 
watershed simulation models, the following criteria may be observed 
for an ideal case: 

1. The model should be sufficiently flexible to allow for an 
interface with other models used in the water management study 
(e.g. optimization of reservoir operations, water quality, 
biomass growth, etc.). 

2. The model structure should be amenable to modifications so as to 
accommodate changes taking place in the watershed, whether they 
be physical or changes in the information being gathered on 
current hydrometeorological conditions. 

3. In addition to the flow simulations, the model should be capable 
of detailed simulations of various other components of the 
hydrologic cycle, facilitating also a total water balance study 
of the basin. 

4. The model should have as part of its simulation package a 
suitable snowmelt routine (for Canadian conditions) so that flow 
simulation during winter months is possible. Alternatively, the 
model should be capable of properly interfacing, on-line, with a 
suitable snowmelt subroutine. 

5. The model should be capable of being calibrated on one or more 
years of daily flows in addition to allowing a calibration on 
seasonal flow data. An ideal case would be a common set or 
grouping of optimized parameters allowing a linking of 
continuous long term and short term simulations. 

6. The model should be programmed to function for different time 
steps, so that short term and long term simulations are possible 
and this should be performed continuously without interruptions. 
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7. The model should be able to simulate flows from the entire 
watershed as well as those from a small sub-basin with a 
reasonable accuracy and without many changes in the parameter 
values. This will facilitate application of the model to 
ungauged sub-basins. 

8. The type of input data (meteorological, physiography of the 
basin, etc.) should be easily obtained by measurements or from 
standard publications. 

9- The model should accommodate both real time simulation and flood 
forecasting requirements. In the latter case, the structure of 
the computer program should be adaptable to modification of 
meteorological data input from other sources such as 
precipitation data from radar systems. 

10. The computer program should be flexible enough to accommodate 
changes in input/output format and to be run on the most readily 
available types of computers. Easy access and use of the 
package program by non-expert hydrological personnel would be an 
advantage. 

The models which are under consideration here do possess some of the 
qualifications mentioned above, but no single model has all these 
properties. The M0EHYDR2 and SQGRID models are partly deterministic 
and parametric and use mathematical relationships based on physical 
concepts in defining various components of the runoff cycle. The 
M0EHYDR2 model gives greater details of various component processes, 
hence the model is useful in water balance studies. The SQGRID 
model gives details of inflow/outflow components and change in 
storage for a particular grid element only. On the other hand, the 
TVA model is based on index methods and does not consider the actual 
physical processes taking place within the runoff cycle. The values 
given by the TVA model for the different components are simply 
indices of moisture allocations by the model. 

The SQGRID model is based on hourly simulations and in combination 
with a large number of grid nodes and extended periods of 
simulation, computer costs would be excessively high. Such a 
distributed model is most suited to study situations where physical 
variations may exist within the watershed and for a case where the 
simulation extends for a few days only. If the model is properly 
interfaced with a precipitation forecasting system, it could perform 
as a useful flood forecasting model. 

The M0EHYDR2 model has the capability for both hourly and daily 
simulation modes with the proper linkage for simulating short term 
and long term flows in a continuous sequence. Because of the nature 
of the model (lumped) in the hourly simulation mode, the prediction 
of peak flows may not be as accurate as those of the SQGRID model. 
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The TVA model uses daily time steps and because of the method of 
calculation (index method), it is useful mainly for studying long 
term simulation. Because of the simplicity of this model, the 
computer costs of its application in simulation for several years of 
flows will likely be quite reasonable and the least of the three 
models considered. 

Because of the above review and qualifications, the models are 
subjectively rated and classified as shown in Table 3«1. In a 
practical situation, however, all of the above good qualities of the 
models are dependent on satisfactory and successful calibration. 
The calibrations of the models are presented in the next chapter. 

3.2 Quantitative Evaluation 

The quantitative evaluations are based on the following criteria and 
considerations. 

(1) Goodness-of-fit statistics: 

- Correlation coefficient (proportion of variance explained by 

calibration) . 

- differences between means and variances 

- residual mass curve 

- comparison of peak flows and times of occurrence, etc. 

(2) Cost of simulation: 

- CPU time for daily simulation 

- CPU time for hourly simulation 
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TABLE 3-1 Qualitative Model Evaluation 



Criteria 



Model Applicability and 
Suitability 



M0EHYDR2 



SQGRID 



TVA 



1. Water balance studies 

2. Water management study 
flow simulation for: 
-days 

-months 
-years 



Good 



Good 
Good 
Fair 



Fair 



Good 



Poor 



Fair 
Good 
Good 



3- Accommodation of local changes 
in watershed characteristics 
-physical 
-hydrometeorological information 

4. Snowmelt routine interface 

5. Continuous long and short 
term simulaton 

6. Transfer to ungauged sub-basins 

7. Input data availability 

8. Flood forecasting 



9. Adaptation to different 
computer systems 

10. Computer costs 



Fair 


Good 


Poor 


Fair 


Fair 


Fair 


Good 


Fair 


Fair 


Good 


Poor 


Poor 


Fair 


Poor 


Fair 


Good 


Good 


Good 



Good 

Good 
Good 



Good Fair 

Good Good 
Poor Good 
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%. MODEL APPLICATION, CALIBRATION AND RESULTS 

4. 1 Test Basins 

The preselected test basins are the Eramosa River, Lutteral Creek 
and Nith River. These basins were selected on the basis of 
representing different physical conditions, in addition to being 
classified as basins with unregulated flows. A map of the Grand 
River basin and the different test basins is shown in Figure 1.1. 

4*1.1 Eramosa River Sub-basin 

The Eramosa River sub-basin has an area of 235 sq. km at the 
location of streamflow gauging station 02GA029 above Guelph. 
Streamflow records are available for the period 1962-1977. Only 
very minor flow regulations are caused by local flow restrictions at 
Everton and Rockwood. The long term mean flow at the gauging 
station is 2.33 m^/s and the recorded maximum and minimum daily 
discharges are 35.1 m3/s in 1972 and 0.1 m^/s in 1964, 
respectively. The length of the Eramosa River sub-basin is about 37 
km and the average width is about 8 km. The river has a length of 
about 40 km with an average grade of 3 m per km. 

4.1.2 Lutteral Creek Sub-basin 

Lutteral Creek is a tributary of the Speed River and the sub-basin 
at the streamflow gauging station 02GA033 near Oustic has a drainage 
area of 65 sq. km. The streamflow is affected by minor local 
restrictions. Continuous flow records of daily discharges are 
available from 1966. The long term mean daily flow is 0.76 nw/s 
and between the period 1966-1976 the recorded maximum and minimum 
daily flows are 13-9 m^/s in 1972 and 0.01 m^/s in 1967, 
respectively. Based on the properties of the Speed and Eramosa 
Rivers, Lutteral Creek may have a grade of about 3 m per km. 

4.1.3 Nith River Sub-basin 

The Nith River sub-basin at the streamflow gauging station 02GA10 
near Canning has a drainage area of 1031 sq. km. The streamflow is 
primarily natural, though because of the meandering nature of the 
stream and gentle slope, a number of naturally occurring oxbows may 
act as local restrictions. Continuous records of mean daily flows 
are available from 1948. A set of previous records are available 
from 1913-1923. The long term mean flow at this location is 10.56 
m3/ s with a maximum daily discharge of 328.5 rn^/s in 1950 and a 
minimum daily discharge of 0.57 m-Vs in 1949 (covering 1948-1976). 

The sub-basin is about 72.4 km in length and 16. 1 km wide. The 
stream has a long meandering course, with an average gradient of 
1.23 m per km for a total length of 158 km. The climatic conditions 
of the sub-basin are very similar to that of the entire Grand River 
basin. 
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GRAND RIVER WATERSHED 
SCALE 

1: 634,000 

EXISTING RESERVOIRS 
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3. CONESTOGA 

4. GUELPH 




▼ Test Basin Gauge Locations 



Figure 4.1 GRAND RIVER WATERSHED AND 
TRIBUTARY AREAS 
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4.1.14 Climate 

There are minor variations in climate throughout the different areas 
of the Grand River basin. The mean annual temperature ranges from 
5.6°C in the north of the watershed, to about 7.8°C in the 
vicinity of Cayuga and the Lake Erie shoreline. The temperature 
during the winter months averages about -6°C with the southerly 
zone slightly warmer than the northern zone. During the summer 
months, the average temperature varies between 17 and 21°C. 

The total precipitation in the watershed during a normal year is 
about 864 mm, out of which the snowfall is about 229 cm or 229 mm 
water equivalent. The seasonal distribution of precipitation tends 
to be regular with the rainfall during May being slightly above the 
average. The northern half of the watershed lies within the 
so-called ' snowbelt' of southern Ontario. Heavy snowfalls in this 
area are common experiences. As a result, the snowpack 
accumulations with accompanying low temperatures, usually persist on 
the ground for long periods. 

4.2 Calibration of Models 

The three models selected are calibrated separately in a single 
sub-basin (Eramosa) with essentially a common set of data. However, 
because of the differences in input requirements for the different 
models, further applications are on a selective basis, depending on 
the availability of common data sets. The M0EHYDR2 and SQGRID 
models are calibrated and compared also in Lutteral Creek sub-basin, 
and the former is calibrated and evaluated further in the Nith River 
basin. 

The calibration involves the determination of an optimal set of 
model parameters. This is achieved by minimizing the sum of the 
squared differences between observed and calculated flows. 



4.2.1 Eramosa River Sub-basin 

The M0EHYDR2 and SQGRID models are calibrated on an hourly basis for 
a common, single storm, rainfall-runoff event while the TVA model is 
calibrated on daily precipitation and runoff data. 

(i) M0EHYDR2 Model Application 

The runoff event selected for the calibration commenced on June 11, 
1971. The same event is used to calibrate the SQGRID model, so that 
a direct comparison of results can be made. 

The records of the rainfall event used cover a period of two days on 
June 13 and 14 with daily total precipitation of 38.6 mm and 19 mm 
for the 1st and 2nd day, respectively. The hourly distribution of 
the rainfall exhibits typically a high degree of variability in 
intensity, with 30.7 mm between 1700 and 1900 hours on the 1st day 
and 14.7 mm between 0300 and 0400 hours on the 2nd day. 
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The basin response to the storm gives a satisfactory recorded output 
hydrograph providing adequate data for the calibration. As only a 
summer rainfall-runoff event is being considered, calibration is 
done for the STSM sub-model parameters. The calibration is 
performed with a combination of coarse trial-and-error runs, with 
refinement using the automatic optimization sub-routine. 

Suitable initial values for the parameters and the initialization 
variables are chosen and the automatic optimization routine is 
allowed to go through 20 iterations. Minor adjustments are made to 
the initial parameter values and the model constants until a 
reasonable fit is obtained between the observed and simulated 
hydrographs. Results of the daily mean of the observed and 
simulated flows and the supporting statistics of fit are given in 
Table 4.1. The optimized set of parameter values are given in Table 
4.2. Suitable estimates of other supporting parameters such as 
those related to ground water storage and flow are determined 
separately by calibration runs on the LTSM sub-model routine. 

(ii) SQGRID Model Application 

The Eramosa sub-basin system is divided into 74 grid elements of 
2 km x 2 km size. The model is then calibrated on the runoff event 
of June 11-18, 1971. In the calibration procedure, suitable initial 
parameter values are chosen and the values for the model constants 
are retrieved from the Basin Files (see Appendix B). Recorded 
hourly flows at the basin outlet are used in the optimizing 
function. The automatic optimization routine "CLIMBX" is allowed to 
go through 16 trials after the initial parameter vaues are adjusted 
until a reasonable fit between the observed and simulateld flows is 
obtained. Table 4.3 shows the optimized parameter values. 

The mean daily values of the observed and simulated flows and the 
statistics of fit are given in Table 4.1 

(iii) TVA Model Application 

The calibration of this model is based on daily precipitation inputs 
and runoff outputs, performed with three years of data. The 
water-years 1968/69, 1969/70 and 1970/71 are used. Calibration is 
done separately for each water-year and finally with the combined 
three years of data. This is done primarily to examine the 
consistency in the values of the estimated parameters with 
different, but similar data sets for a given basin. Table 4.4 shows 
a summary of the parameter values obtained for the various 
calibrations. It is evident that there is a significant range for 
each parameter. The lack of consistency in parameter values 
indicates the optimized set of parameters are highly unique for each 
data set and may not be satisfactory for long term simulation. 
Further refinement in the calibration may result in improvement. 

Table 4.5 gives the comparison of observed runoff with the simulated 
values for the three years 68/69, 69/70 and 70/71. The results 
indicate that total yearly runoff is simulated reasonably well. 
However, during low flows (monthly) the model simulates flows higher 
than the observed; the reverse holds for the high flow simulation 
periods. 
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TABLE 4.1 



Observed and Simulated Flows Using Models M0EHYDR2 and SQGRID 
Eramosa River Sub-basin 



Date 


Observed Flow Simulated (Calibrated) Flow 




M0EHYDR2 (20 iterations) SQGRID (16 iterations) 


June '71 


m3/s m^/s m3/s 



11. 


1.21 


12. 


1.03 


13. 


1.13 


14. 


4.28 


15. 


9.20 


16. 


6.97 


17. 


3.23 


18. 


1.70 



Mean 3- 59 

Std. Dev. 3.06 

Correlation 
Coefficient 

CPU time 



1.19 


0.02 


1.10 


0.07 


1.15 


0.99 


6.92 


4.39 


7.43 


7.26 


4.77 


8.99 


3.16 


5.72 


2.19 


4.37 


3.49 


3.98 


2.60 


3.36 


0.876 


0.848 


24.9 sec 


2 min 17 sec 



TABLE 4.2 



Optimized Parameter Values for the M0EHYDR2 Model 
Calibration: Eramosa River Sub-basin (Event: June 11-18, 1971) 



Parameter* 



Optimized Value 



A 

AN 

AL 

NN 

FINC 



10.99 mm/hr 
3.82 
1.05 
0.19 
4.21 mm/hr 



* See Appendix A for description 
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TABLE 4.3 



Optimized Parameter Values for the SQGRID Model Calibration : 
Eramosa River Sub-basin (Event: June 11-18, 197 1 ) 



Parameter* Optimized Value 





FCFRAC 




3.63 




GWFRAC 




0.002 




SFCLAG 




2.4 




SFOLAG 




1.71 




* See Appendix 


B 




TABLE 4.4 









Optimized Parameter Values for the TVA Model Calibration : 
Eramosa River Sub-basin (Event: 1968-1971) 



Parameter* 
Water Year A D GWK TDSRO 



1968-1969 


12.23 


96.1 


5.75 


.182 


.751 


1969-1970 


56.00 


26.0 


5.70 


.075 


.327 


1970-1971 


4.53 


74.19 


14.25 


.168 


.272 


1968-1971 


3.87 


52.50 


13-47 


.153 


.286 



* See Appendix C 
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TABLE 4.5 



Observed and Simulated Runoff-TVA Model-Daily Simulation : 











Eramosa River Sub 


-basin. 










68-69 




69-70 




70-71 




Month 


Observed 
mm 




Simulated 
mm 


Observed 
mm 


Simulated 
mm 


Observed 
mm 


Simulated 


October 


19.23 




35.31 


7.32 


12.75 


16.76 


17.65 


November 


30.48 




56.39 


19.20 


36.07 


24.16 


37.53 


December 


31.24 




47.75 


11.33 


19.05 


34.80 


43.95 


January 


29.97 




24.71 


9.35 


2.21 


18.64 


4.52 


February 


30.73 




18.44 


9.93 


1.17 


19.99 


1.42 


March 


62.48 




53.85 


15.44 


7.14 


32.00 


23.93 


April 


92.46 




72.39 


71.37 


52.07 


87.38 


71.38 


May 


43-94 




16.94 


28.19 


21.41 


22.61 


21.87 


June 


16.79 




12.90 


10.59 


16.43 


22.22 


26.67 


July 


9-93 




12.12 


8.81 


14.86 


16.76 


27.69 


August 


8.99 




14.15 


9.68 


14.83 


18.39 


30.73 


September 


4.47 




10.69 


13.31 


12.90 


11.68 


25.40 


Total 


380.75 




375.67 


214.63 


211 .33 


325.12 


333.25 



4.2.2 Lutteral Creek Sub-basin 

(i) M0EHYDR2 Model Application 

The hourly simulation mode of M0EHYDR2 model is calibrated on the 
runoff event of October 20-27, 1972. The observed rainfall on the 
2nd, 3rd and 4 th days amounts to 8.4, 28.5 and 31.8 mm, 
respectively. The distribution of the rainfall intensity was 
approximately constant in the order of 2.5 mm/hr during the 3rd and 
the 4th days. No snowfall was recorded during the event and the 
potential evaporation varied from 1.02 mm to 4.6 mm/day. The 
optimization routine is specified with a maximum of 20 iterations 
and Table 4.6 gives the optimized parameter values. 

Table 4.7 gives the daily mean values and the summary statistics for 
the observed and simulated flows. 

(ii) SQGRID Model Application 

The Lutteral Creek sub-basin is divided into 16 grid elements of 
2 km x 2 km size and the model is calibrated on the runoff event of 
October 20-27, 1972. Table 4.8 gives the initial and the optimized 
parameter values. The mean daily values of the observed and 
simulated flows and the statistics of fit are given in Table 4.7. 
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TABLE 4.6 



Optimized Parameter Values for the M0EHYDR2 Model Calibration : 
Lutteral Creek Sub-basin (Event: October 20-27 , 1 972 ) 



Parameter* 



A 

AN 

AL 

NN 

FINC 



Optimized Value 
2.58 mm/ hr 
2.38 
2.60 
0.27 
0.397 mm/hr 



* See Appendix A 



TABLE 4.7 



Observed and Simulated Flows Using Models M0EHYDR2 and SQGRID 

Lutteral Creek Sub-basin 



Date 


Observed Flow 


Simulated (Calibrated) Flow 




■3/a 


M0EHYDR2 (20 iterations) SQGRID (16 iterations) 


Oct. 72 


m3/s m3/s 



20 


0.12 


21 


0.13 


22 


0.27 


23 


1.77 


24 


2.05 


25 


1.42 


26 


0.99 


Mean 


0.92 


St.Dev. 


0.76 


Correlation 




Coeff. 




CPU time 





0.12 
0.13 
0.15 
0.88 

2.39 
1.50 
0.92 

0.89 
0.85 

0.87 

22.5 sec 



0.08 
0.26 
0.50 
1.51 
2.07 
1.15 
0.97 

0.94 
0.66 

0.96 

2 min, 
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4.2.3 Nith River Sub-basin 

M0EHYDR2 Model Application 

The M0EHYDR2 model is calibrated on daily streamflow data for the 
period October 1 - December 31, 1 974 . This calibration is done for 
the daily simulation mode (LTSM Sub-model). In this mode the 
parameters which are sensitive to flows, are different from those in 
the hourly simulation mode. Seven model parameters are optimized 
after establishing suitable initial values by preliminary and 
supplementary hydrologic analyses. Refinements in these values are 
achieved with the automatic optimization sub-routine. 

The precipitation data for this period indicated rainfall of 42.2, 
87.1, and 8.9 mm, respectively, occurred during October, November 
and December. The corresponding snowfall amounts are 3.3, 20.3 and 
42.4 mm. The average temperature during October was above 0°C and 
that during November was about -0.5°C. Temperatures during 
December were always below 0°C. The optimized parameter values 
are given in Table 4.9- The corresponding values of the model 
constants are given in Appendix A. 

Figure 4.2 shows the plot of observed vs. simulated (calibrated) 
flows for the period under consideration. A visual examination 
shows that the simulation is reasonably good. The mean and standard 
deviation of the observed data are 4.13 m3/s and 1.95 m^/s 
respectively and that of the simulated/calibrated data are 4.33 m^/s 
and 1.62 m3/ s . The coefficient of correlation is 0.82. The 
scatter diagram as presented in Figure 4.3, shows that there is a 
tendency for underestimation of high flows, and a tendency for 
overestimation of low flows, though less significantly. Further 
refinements in the model parameters, primarily those dealing with 
baseflow calculations, through a larger number of optimization runs 
and a better choice of initial values could improve the fit. 

Figure 4.4 shows a plot of the residual mass curves of the observed 
and simulated flows. The residual mass curve is the cumulative 
difference between the individial flow (daily) and the mean flow. 
Thus, the ordinate of the residual mass curve, at any point in time, 
depends on the history of the preceeding events. Any deviation of 
the residual mass curve of the estimated flows from that of the 
observed flows would indicate the existence of systematic errors in 
the estimated flows. The plots show a reasonable fit during the 
first half of the simulation period and that systematic errors are 
introduced during the second half period. Systematic errors are 
common in rainfall-runoff hydrologic models and a simple test to 
detect these errors is important in water resources engineering 
(Aitken, 1973). 
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TABLE l*. 8 



Optimized Parameter Values for the SQGRID Model Calibration : 
Lutteral Creek Sub-basin (Event: October 20-27, 1972) 



Parameter* Optimized Value 

FCFRAC 1.027 

GWFRAC 0.0138 

SFCLAG 0.919 

SFOLAG 2.425 



* See Appendix B 



TABLE 4.9 



Optimized Parameter Values for the M0EHYDR2 Model 
Daily Simulation - Nith River Sub-basin 
(Event: Oct.1 - Dec. 31, 197 l O" 



Parameter* Optimized Value 

CP 71.12 mm/ day 

FINC/SHC 24.38 mm/ day 

EP 0.45 

EM 0.40 

C1 0.29 

CS 0.049 

TRANS 0.004 mm/hr 



* See Appendix A 
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5. MODEL PREDICTION RESULTS 

In Chapter 4 the optimized parameter values for the three models, as 
applicable to the different basins, are given. In this Chapter, the 
optimized parameters are used to predict flows for a different 
period. 

5.1 M0EHYDR2 Model Predictions 

The optimized parameter values for the Nith River sub-basin, as 
given in Table 4.9, are used to predict daily flows for the period 
May-July, 1974. The meteorological data for this period include 
rainfall totals of 124.7 mm, 86.4 mm and 58.7 mm for the months of 
May, June and July, respectively. There is no record of snowfall 
during this period. The peak discharge on May 18th results mainly 
from a rainfall input of 47.5 mm on May 16th. The predicted flows 
are compared with the observed flows in Figure 5.1. The predicted 
hydrograh retains the general shape of the observed hydrograph, 
though the peak flow is underestimated (and hence underestimates the 
mean flow). There appears to be less than satisfactory predictions 
of observed minor fluctuations during the low flow period. However, 
these flow fluctuations could have been accentuated mainly by 
artificial and local channel storage effects rather than true basin 
response to precipitation and potential evapotranspiration input. 
Futhermore, a single period model calibration is not enough to 
produce sufficient statistics of the model parameters for use in 
predictions. 

5.2 SQGRID Model Predictions 

The optimized parameter values for the Eramosa sub-basin in Table 

4.3 are used to predict flows for the runoff event of August 25 - 
September 1, 1971. The mean daily observed and predicted flows are 
compared in Table 5.1. It can be seen that the predicted low flows 
do not agree very well with the observed flows. The predicted high 
flows agree well with the observed values, but lag them by about one 
day. A shift in peak is also observed in the calibration process. 
It is therefore presumed that further refinement in calibration may 
improve the timing of the predicted peak flows. 

5.3 TVA Model Predictions 

The TVA model, with the optimized parameter values of Table 4.4, is 
used to predict flows in the Eramosa River and the Speed River. The 
Eramosa River is a tributary of the Speed River and it is assumed 
that the Speed River basin has similar physical characteristics as 
those of the Eramosa River sub-basin. The optimal set of parameter 
values for the calibration based on the combined three years of data 
is used (Table 4.4). The observed and predicted monthly runoffs in 
the Eramosa River for the water year 1971-1972 are given in Table 
5.2. The observed and predicted runoffs exhibit relations similar 
to those indicated by the results of Table 4.5. The yearly total 
runoff is predicted with reasonable accuracy. 
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TABLE 5.1 

Observed and Predicted Flows-SQGRID Model 
Eramosa Sub-basin (Event: August 25 - September 1, 1971 ) 



Predicted Flow 



Date 






Observed Flow 


Aug-Sept 


•71 


m3/s 


25 






0.71 


26 






0.89 


27 






6.67 


28 






8.77 


29 






7. 44 


30 






4.40 


31 






2.88 


01 






1.71 


Mean 






4.12 


Std. 


Deviation 


3.13 


Corr. 


coe 


fficient 




CPU 1 


'ime 







m3/j 



0.565 



0.05 
0.13 
1.72 
3-94 
8.34 
6.61 
4.49 
3.23 

3.56 
2.95 

8.05 Sec 



TABLE 5.2 

Observed and Predicted Runoff-TVA Model 
Eramosa Sub-basin (October 1971 - September 1972 ) 

Month Observed Predicted 



mm/month mm/month 



October, 71 9.19 18.11 

November, 71 11.66 25.10 

December, 71 27.94 48.26 

January, 72 18.80 15.62 

February, 72 15.80 3-66 

March, 72 24.84 45.97 

April, 72 132.08 60.45 

May, 72 37.34 35.05 

June, 72 25-65 34.29 

July, 72 19.41 29.97 

August, 72 8.97 22.02 

September, 72 8.66 16.84 



Mean 28.45 29.46 

Std. Dev. 33-78 16.13 

Corr. Coeff. 0.719 

CPU Time 0.2 Sec. 
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In the Speed River basin, the monthly values of the observed and 
predicted runoffs for the water years 1969-1970, 1970-1971, and 
197 1-1 972 are plotted in Figure 5.2. A visual inspection of the 
graphs shows that the predicted runoff compares reasonably well with 
the observed values in many months. Table 5.3 summarizes the 
statistical analysis for examining the significance of the 
difference between the observed and predicted flows at the 95 per 
cent confidence level. The results indicate that for the months 
February, May, September and October, the predicted flows are 
statistically significantly different from the observed at the 
specified significance level. 

Based on the limited model prediction results reported herein, the 
qualities of prediction by the three models are summarized in Table 
5.^., together with information on CPU costs. 
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TABLE 5.3 



Statistics of Observed and Predicted Runoff-Speed River 
1970-1971 Flows - TVA Model Prediction 



Observed 



Predicted 



Mean 
mm/day 



Std. Dev, 
mm/day 



Mean Std. Dev. Computed 
mm/day mm/day t 



October 


0.515 


0.203 


0.358 


0.152 


3-45 


November 


0.808 


0.305 


0.846 


0.178 


0.59 


December 


1.054 


0.533 


1.303 


0.686 


1.59 


January 


0.589 


0.178 


0.467 


0.406 


1.53 


February 


0.617 


0.127 


0.246 


0.432 


4.35 


March 


1.069 


0.330 


1.405 


1.346 


1.35 


April 


3.132 


1.702 


2.964 


2.362 


0.315 


May 


0.650 


0.305 


0.391 


0.051 


4.67 


June 


0.724 


0.660 


0.543 


0.279 


1.38 


July 


0.531 


0.432 


0.63 


0.178 


1.18 


August 


0.592 


0.787 


0.671 


0.152 


0.56 


September 


0.391 


0.152 


0.577 


0.051 


6.42 


TABLE 5.4 
















Quality of Model 


Predict: 


Ions 





Criteria 



Model 



M0EHYDR2 



SQGRID 



TVA 



1. 


Mean monthly 


flow 


Good 


— 


Good 


2. 


Mean daily fl 


ow 


Good 


Good 


Fair 


h 


Peak flows 












- hourly (storms) 


Good 


Good 


— 




- daily 




Fair 


— 


Fair 


n. 


Low flows 




Good 


Poor 


Fair 


5. 


CPU costs for 
predictions 




Moderate 


Highest 


Lowest 
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6. DISCUSSION OF RESULTS 

The model calibrations presented in Chapter 4, are carried out on 
one data set for each model. Tests for model predictions are 
demonstrated in Chapter 5. In a practical situation, the models 
should be calibrated on data sets from different years. The average 
values for the parameters determined in these calibrations should be 
used in predictions. The several calibrations are necessary because 
the optimized parameter values vary significantly with the different 
sets of initial parameter values chosen, the model constants and the 
types of data sets used. These variations in optimized parameter 
values are particularly significant in complex models where indices 
of many hydrologic processes are interrelated with one another. The 
above phenomenon can be observed in Table 4.4 in which the sets of 
optimized parameter values obtained from the calibration of the TVA 
model are quite different from one set to another. However, in 
Table 4.5 all three cases of simulation give reasonable results with 
respect to the annual mean runoff. 

The results of model calibration for the Eramosa sub-basin presented 
in Table 4.1 indicate that the M0EHYDR2 simulated mean flow and the 
standard deviation of flows are within three per cent and 15 per 
cent of the observed flow values, respectively. The corresponding 
flow statistics values for SQGRID simulations are 11 per cent and 10 
per cent respectively. The coefficients of correlation between 
observed and simulated flows are 0.876 for M0EHYDR2 and 0.848 for 
SQGRID, respectively. This shows that the simulations are generally 
good. However, the M0EHYDR2 model underestimates the peak by about 
19 per cent and the SQGRID peak flow lags the observed peak by one 
day. Furthermore, the poor simulations by SQGRID during the initial 
period may be attributed to the improper accounting of the initial 
baseflows. The CPU time of 24.9 sec (20 iterations) for the 
M0EHYDR2 simulation as compared to 2 min 16.9 sec (16 iterations) 
for the SQGRID simulation makes the M0EHYDR2 model more attractive, 
depending upon the required accuracy for simulation. 

The shortcomings observed above in the results of simulation of 
Eramosa River flows do not occur in the results of simulation of the 
Lutteral Creek flows. In Table 4.7, it can be seen that the peak of 
SQGRID simulated flows coincides with the observed peak and the 
M0EHYDR2 simulated peak is 16.6 per cent higher than the observed. 
The coefficient of correlation between observed and simulated flows 
remains constant for M0EHYDR2, whereas that of SQGRID increases to 
0.96. The above calibration shows that the errors introduced may 
not be a function of the modelling technique itself. Improved 
simulations can be obtained by a proper choice of initial values for 
parameters and model constants. Further efforts with the automatic 
optimization may improve the result of calibration. 

The results of simulation of Nith River flows by the daily 
simulation mode of M0EHYDR2 presented in Figure 4.1, is more 
encouraging than the results obtained by the hourly simulation mode 
of M0EHYDR2, when applied to Eramosa River and Lutteral Creek 
sub-basins, as presented in Tables 4.1 and 4.7. The point that 
should be noted here is that, with hourly flows, the accuracy of 
simulation should be greater to get comparative results with that of 
daily simulations. Furthermore, in the hourly simulated flows, the 
number of data points available for optimization (8 points in cases 

here) is less than that for the daily simulation (92 points for Nith 
River) . 
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The predicted daily flows in the Nith River, presented in Figure 
5.1, are reasonably good, noting that the model parameters are 
obtained after calibration on one data set. The reason for the 
underestimation of the peak and mean flows may be due to excessive 
evaporation rates during the peak flow period. The evaporation 
rates are based on records from the Elora research station which is 
completely outside the Nith River sub-basin. A check on the water 
balance shows that for the month of June, about 84 per cent of the 
total rainfall is lost to evaporation. 

An attempt is made also to predict the Nith River flows for one full 
year, based on the optimized parameter values. However, the 
predictions (not presented here) are far from satisfactory. This 
indicates that prediction for longer time periods should be based on 
calibrations which have data sets of the same time period or longer. 

An attempt to calibrate the M0EHYDR2 model for the Nith River 
sub-basin on one full year's data set (October 74 to September 75) 
was not successful. The calibration procedure encountered 
difficulties in the selection of a good set of initial parameter 
values, so that the optimized values would give flows which are 
comparable to the observed flows. The drainage area of the Nith 
River basin, 1031 sq. km, is approximately a fifth of the entire 
Grand River watershed. The parameters which undergo optimization in 
the M0EHYDR2 model each represent some physical characteristic of 
the basin, and may have significant variations spatially and 
temporally. The selection of a set of parameter values to represent 
the average conditions of the basin, for all seasons of the year, 
therefore, becomes an important and very difficult task. After 
several trial simulations, it is found that the simulated flows 
compare very well with the observed flows only during a portion of 
the year (May-September). Since the parameter values are calibrated 
on one full year's observed flows, they tend to represent average 
conditions of the basin resulting in an underestimation of peak 
flows and overestimation of low flows. The means of the observed 
and simulated flows are, respectively, 11.8 and 9.5 m^/s which 
shows the overall mean is underestimated. However, if March and 
April (months which have high observed peaks, but are severely 
underestimated) are omitted, then the corresponding values are 10.2 
and 10.6 m3/s, which shows an overestimation. This indicates that 
if March and April are omitted, a better simulation may be obtained. 

These calibrations show that in large basins, the M0EHYDR2 model is 
easily calibrated when the simulation period is seasonal. 
Furthermore, it is found that the effort necessary (trial runs, 
time, etc.) for calibration, is greater when the data set used in 
calibration has a longer period. 

On a yearly data set, the effort necessary to calibrate the TVA 
model is less than that in case of the M0EHYDR2 model. This is 
because the parameters of the TVA model are based on empirical 
equations and they do not represent any actual characteristics of 
the basin. However, in the TVA model, this may lead to incorrect 
accounting of the different components (eg. baseflow) of the runoff 
cycle, even though the total runoff is correctly simulated. 
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The results presented in Table 5.2, the summary of SQGRID model 
predicted hourly flows in the Eramosa River, are comparable to the 
results of calibration given in Table 4.1. The predicted peak is 
within 5 per cent of the observed peak even though the former one is 
lagging by one day. This lagging of daily peak flow is also noted 
in the calibrations. The poor simulations at the beginning of the 
period may be due to incorrect baseflows which are not considered in 
the model. 

The comparison of observed runoff with TVA model predictions in the 
Eramosa River, presented in Table 5.3, shows that there is a good 
agreement between the annual mean flows. However, there are 
significant differences in the monthly mean values for four months. 
Imprecisions are reflected also in some of the standard deviations. 
Based on a daily basis the discrepancies may be higher. The results 
of predicted runoff in the Speed River, plotted in Figure 5.1, show 
that the difference between observed and predicted runoff are of the 
same order of magnitude as that in the Eramosa sub-basin. The poor 
simulations of the TVA model on a daily basis are revealed in the 
statistical analysis (Student's t test) given in Table 5.3. This 
evaluation shows the need for statistical analysis of the observed 
and predicted flows in order to assess the quality of simulation by 
a model. 
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7. CONCLUSIONS 

The three watershed models M0EHYDR2, SQGRID and TVA are evaluated 
both qualitatively and quantitatively for their applicability in the 
Grand River basin. The qualitative evaluation is given in Chapter 3 
whereas Chapters 4, 5 and 6 deal with the quantitative evaluation 
aspects. 

i A direct comparison between the SQGRID model and the M0EHYDR2 

model for hourly simulations shows that the models yield results 
which appear to have the same degree of accuracy. The M0EHYDR2 
model is more attractive because of less CPU time involved. In 
Chapter 3, it is observed that the SRGRID model could be used as 
a good flood forecasting model. However, in the case of the 
Grand River basin, the history of flooding reveals that the 
flooding is mainly due to the snow melting associated with 
spring rain. Therefore, the SQGRID model can be used to 
forecast springflood only by interfacing the model with a proper 
snowmelt routine. Since the ground surface and the top soil 
conditions during spring will be different from conditions in 
any other season, it is essential that the model, interfaced 
with a snowmelt routine, be calibrated on several spring storm 
flooding events if it is to be used for forecasting. 

The M0EHYDR2 model is recommended for daily flow simulations. 
The SQGRID model cannot be used as a daily simulation model 
because the time step is restricted by stability criterion based 
on the grid size. However, the model is useful to study the 
effect of changes in the distribution of certain specified 
physical characteristics of the watershed. 

ii The TVA model, the least complex model, has small computer time 
requirements. It is recommended for long term simulations where 
average monthly flows will satisfy the flow requirements. The 
TVA model is the only model in which the optimized parameter 
values can be obtained based on calibration using continuous 
flow data over several years, and hence will give parameters 
which are temporally averaged. 

iii The M0EHYDR2 model is the one model which has the hourly as well 
as the daily simulation modes, and hence is applicable in a 
variety of situations in a comprehensive water resources 
management program. Furthermore, a snowmelt subroutine is a 
component of the M0EHYDR2 model. 
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APPENDIX A 
M0EHYDR2 MODEL 
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APPENDIX A - M0EHYDR2 MODEL 
(i) LTSM Sub-Model - Parameters Optimized 

CP = A multiplication constant used in calculating 

infiltration 
EP = An exponent for precipitation used in calculating 

infiltration 
EM = An exponent for soil moisture used in calculating 

infiltration 
TRANS = Transmissibility factor for baseflow calculations 
CS = Constant for determining change in depth of water in 

baseflow storage (index of storage coefficient) 
FINC = Steady state infiltration rate 

SHC = Steady state percolation rate from interflow zone 
C-i = Critical ratio of volume of water in upper half of 

layer to storage in total interflow zone layer 



(ii) 



FINC 



AN 
AL 



NN 
(iii) 



STSM Sub-Model - Parameters Optimized 

= Steady state infiltration rate 

= A multiplication coefficient used in calculation of 

infiltration 
= An exponent used in calculation of infiltration 
= A multiplication coefficient used in determining 

surface runoff 
= An exponent used in determining surface runoff 

Model Constants 



AREA = Drainage area of basin 

ISCM = Maximum interception storage capacity 

FARIS = Fraction of area assigned to interception storage 

DSCM = Maximum depression storage capacity 

FARDS = Fraction of area assigned to depression storage 

TS = Total storage in the profile 

NLAYER= Number of layers in profile not including interflow 

zone 
FCAP = A vector containing the field capacities of each layer 

in the profile 
WP = A vector containing the wilting points of the layers 

of the profile 
ALARGL= Length of stream 

ALITL = Average distance of drainage basin to center of stream 
HMIN = Minimum depth of water in baseflow storage 
K1 = Routing constant for overland flow 
K2 = Routing constant for interflow 
K3 = Routing constant for baseflow 
B = Coefficient used in calculation of evapotranspiration 

from soil profile 
RSRO = Value of routed surface runoff rate 
HACT = Depth of water in baseflow storage 
M = The number of points in the time-area histogram for 

lagging channel flow in STSM 
TARHIS= Time-area histogram for lagging channel flows in STSM 
AAP = Regression constant in snow density-depth relation 

(accumulation period, no rain) 
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BAP = Regression coefficient in snow-density-depth relation 

(accumulation period, no rain) 
ARAP = Regression constant in snow density-depth relation 

(accumulation period, rain) 
BRAP = Regression coefficient in snow density-depth relation 

(accumulation period, rain) 
AMP = Regression constant in snow density-depth relation 

(melt period, no rain) 
BMP = Regression coefficient in snow density-depth relation 

(melt period, no rain) 
ARMP = Regression constant in snow density-depth relation 

(melt period, rain) 
BRMP = Regression coefficient in snow density-depth relation 

(melt period, rain) 
BRAP1 = Regression coefficient in snow density-depth relation 

(accumulation period, rain) 
BRMP1 = Regression coefficient in snow density-depth relation 

(melt period, rain) 
BS = Constant used in calculating fraction of interception 

storage in SNMELT 



(iv) 



Printout Options 



Option 1 Table of monthly summary for the different hydrologic 

component blocks. 
Option 2 Table of simulated daily flows. 
Option 3 Table of output for the snowmelt sub-model. 
Option 4 Table of daily summaries of infiltration and soil 

moisture storages. 
Option 5 Table of daily output and storages for interflow and 

the baseflow components. 
Option 6 Table of daily summaries of the primary hydrologic 

components of the simulation. 
Option 7 Plots of precipitation, soil moisture storages, actual 

evaporation, observed flows and simulated flows. 
Option 8 Plots of simulated snowpack depth, water equivalent 

and density. 
Option 9 Summary of hourly simulations (outputs and storages). 
Option 10 Summary of model parameters (including optimized 

values) . 
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month; 11 



UAY 



10 
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14 
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0.0 
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_0*JL 
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*X 
JUL 



*X 
JUL 
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0.0 1 / fl.no/ 0.01 / «.oi/ 



* SM«fcLl - S^»i*MtLI «Ot»El 
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(i) 



APPENDIX B - SQGRID Model 

Structure of the Computer Program 
MAIN 



ERRSET CLIMBX 



DRBRMV 



fcn i r v r n r I | 

OLAG RETVE GETMET INFILT SETMET DAYEVP CLAGNS ; 

i CLAGAS 

EVTURC 

MAIN = Main program 

ERRSET= Statistics of errors between observed and simulated 

flows 
DRBRMV= Runoff model 
FCN = Optimizing function 

GETMET= Extract (reformat) meteorological data 
CLIMBX= "hill-climbing" optimization routine 
OLAG = Computes overland lag (in hours) 
RETVE = Reads records from basin file 
INFILT= Computes infiltration from surface into soil, and 

percolation from soil to ground water storage 
SETMET= Reads temperature and precipitation data 
DAYEVP= Calculates daily evaporation 

EVTURC= Calculates annual evaporation using Turc's formula 
CLAGNS= Computes channel lag (in hours) for natural conditions 
CLAGAS= Computes channel lag (in hours) for artificial 

conditions 



(ii) 



Basin File (BASFI) 



The Basin File (BASFI) is a direct access file containing 33 records 
of 2730 half word integers for the 2730, 2 x 2 km grid cells in the 
study area. Each record represents a particular physical 
characteristic of the study area. Data in the file are stored by 
column, from bottom to top and from left to right. The Grand River 
above Paris is represented by a 39 x 70 matrix, or 2730 cells. 



70 



■t-4- 



39 
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Layout of the Fil e 

ftfefcftrd Content? 

riv. ie) 

t ''low (1=nortn, 2ceast, 3 2 scuth, 

'J=west, Osnot in atudy area) 

4 Number of squares water has flowed through to 
this square 

5 Percentage of infiltration capacity: 

total # land uses 

j (L^Z land use i)*100 

where Lj_ represents the infiltration capacity 
of land i 

6 Maximum voids in soil (%) 

# of soil types 

/ , (Sj*$ soil j)*100 

j=1 in grid 

(80* of this figure is the maximum soil moisture 
available for evaporation) 

7 Constant rate of infiltration after prolonged 
wetting (inches/hour * 1000) 

8 Percentage of soil Type 1 - till moraine 

9 Percentage of soil Type 2 - till plain 

10 Percentage of soil Type 3 - clay plain 

11 Percentage of soil Type 4 - sand plain 

12 Percentage of soil Type 5 - limestone plain 

13 Percentage of soil Type 6 - swamp and bog 

14 Percentage of soil Type 7 - karae moraine 

15 Percentage of soil Type 8 - water 

16 Percentage of soil Type 9 - Canadian Shield 

17 Percentage of soil Type 10- Niagara Escarpment 

18 Percentage of land use Type 1 - industrial 

19 Percentage of land use Type 2 - residential 

20 Percentage of land use Type 3 - associated urban 

21 Percentage of land use Type 4 - horticulture 

22 Percentage of land use Type 5 - orchard and 
vinyards 

23 Percentage of laud use Type 6 - hay, grain and 
soybeans 

24 Percentage of land use Type 7 - tobacco 

25 Percentage of land use Type 8 - improved pasture 

26 Percentage of land use Type 9 - open and scrub 
grassland 
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Record 



Contents 



27 
28 
29 
30 

31 

32 
33 



Percentage of 
Percentage of 
Percentage of 
Percentage of 
Percentage of 
over woodland 
Percentage of 
Percentage of 



land 


use 


Type 


10 


land 


use 


Type 


11 


land 


use 


Type 


12 


land 


use 


Type 


13 


land 


use 


Type 


14 


land 


use 


Type 


15 


land 


use 


Type 


16 



unproductive land 
dense woodland 
open woodland 
scrub woodland 
cut over burnt 

swamp and marsh 
lakes and rivers 



(iii) Program BASFIA 
Purpose 

Creation of four records of a direct access basin file, 
determination of the hierarchial order and establishment of the 
drainage basin, using "manually defined" flow directions. 

Preparation of Input for BASFIA 

(a) Choice of suitable grid size: 

The grid size refers to the length of one side of a grid element on 
the earth's surface. There are several criteria to determine the 
size of the grid. 



1. 



Level of detail available in information (i.e. the scale of 
topographic maps, soil maps, land use maps, etc.). 



General Guide: 

Scale 

1:50,000 
1:25,000 



Grid Size 

5-10 km 
1-2 km 



2. 



3- 



Required level of detail in simulation - The purpose of the 
study must be considered. A smaller grid size gives a higher 
degree of accuracy. 

Financial constraints - The time, budget, and capability of the 
computer system are important. Decreasing the grid size by half 
could increase work load, CPU time, core requirements and disk 
storeage by four to six times. 

(b) Establish grid system and direction of flow: 

On a small scale topographic map, the boundary of the river basin 
and the direction of flow for each grid is usally well defined by 
the channel network and elevation. 
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Flows are coded as follows: 



(iv) 



Purpose 



h 



not in study area 
flow to north 
flow to east 
flow to south 
flow to west 



Program BASFIB 



To overlay digitized spatially distributed data on soil types and 
land use to the basin file. 



Preparation of Input : 

Digitizing and editing - digitize the land use/ land cover map or the 
soil map covering the study area. Edit the digitized data and 
create a VBS file to store the edited data as follows: 

ICODE, NPTS, NPOL, (X(I), Y(I), 1=1, NPTS) 

where ICODE= 



NPTS = 



code of land use or soil (i.e. the record 
number in the basin file) 

number of digitized points in particular 
polygon 



NPOL = sequence number of the polygon 
X(I), Y(I) = Co-ordinates of the vertex of the polygon 
Program BASFIC 



Purpose : 



Delineation of a sub-basin above any given point on the 39 x 70 
matrix which defines the Grand River study area. 



Input: 



Unit 5 
ICOUT = 



column number of the "mouth" of the 
sub-basin being considered (a number from 
1-39) 



JROUT = row number of the "mouth" of the sub-basin ( 
a number from 1-70) 

(vi) Program BASFID 

Purpose : 

To generate both digital and gray maps showing the spatial 
distribution of data in the basin file (eg. soil, land use data) 

Input : 

One card is input for each set of maps required. The card adopts 
the following format: 

ISOIL, TITLE (format: 12, 8A8) 

ISOIL = record number which is used to store desired 
spatially distributed data (e.g. 8 = till moraine, 
21 = horticultural land use) 

TITLE = title or name of record 

(vii) Program BASFIE 

Purpose : 

To calculate the following information, using soil, land use, and 
land-cover distribution. 

a - principal land use (record 5 of basin file) 

b - maximum soil moisture available for evaporation (record 6 

of basin file):80% of maximum voids in soil = maximum soil 

moisture available for evaporation 
c - constant rate of infiltration after prolonged wetting 

(record 7 of basin file) 

Input : 

Input from basin file only: no card image input required. 

(viii) Input Parameters and Variables 

Card No. Format Contents 

1 20A4 Title 

2 10F5.0 SAINIT, FCFRAC, GWFRAC, 

GWINIT, GWLAG, SFOLAG, 
SOILT, SNOWT, TIMEIN 

3 815 NSQ, INYR, INMTH, INDAY 

INDAYS, MASDIV, IOPT, 
NPRIME 

4 215 NIRIPV, NVAR 
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5 

6 

7 

8-n 

n-end 



4F10.0 

4F10.0 

215 

214 

6F10.0 



G(I), 1=1, NVAR 
H(I), 1=1, NVAR 
NLOOK, LOOKSQ 
NREC, NXISQ 
QMESD 



SAINIT = Initial available soil moisture storage 

FCFRAC = Infiltration scaling factor 

GWFRAC = Fraction of ground water flow to base flow 

SFCLAG = Travel time of channel flow scaling factor 

GWINIT = Initial ground water storage 

GLWAG = Travel time for ground water flow 

SFOLAG = Overland flow travel time scaling factor 

SOILT = Soil thickness 

SNOWT = Snow thickness 

TIMEIN = Basic time interval 

NSQ = Number of squares in sub-basin 

INYR,INMTH,INDAY 

= beginning year, month, day of simulation 

INDAYS = Number of days in simulation 

MAXDIV = Subdivisions of TIMEIN 

IOPT = Set to 1 for CLIMBX optimization 

NPRIME = 999 for brief output 

NIRIPV = Number of trials in CLIMBX 

NVAR = Number of model parameters being optimized 

G = Lower bounds of model parameters being optimized 

H = Upper bounds of model parameters being optimized 

NLOOK = Number of squares for which detailed output is 

required 

LOOKSQ = Subscript number of the squares for which 

detailed output is required 

NREC = Location in basin file record of each square 

NXISQ = Number of square water flows to from each square 

QMESD = Hourly flow data 
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ITR DAY MR OV PREC 

IN. 



TE^P EVAP SFSTOR SFlNCl OVLAC CMLAC 3F0UTA W SAl OTFC GhlNCl G*5T0R GwOyTA 
D FAR IN, IN. IN, MR, MR, IN/MR IN, IN, IN, IN, IN, Im/h« 



en 



145 
146 
14T 
140 
149 
150 
151 
152 
153 
isa 

155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 



1 

2 

3 

a 

5 

6 

7 

R 

9 

10 

11 

12 

13 

14 

15 

16 

17 

IP 

1* 

20 

?1 

22 

23 

21 



1 0,0 

1 0.0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 o t o 

1 0.0 

1 0.0 

1 0.0 

1 0.0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 



clock n 



7) « 



169 

170 
171 

.17* 
173 

174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 



8 
8 
8 
8 
8 
8 
8 
8 
8 

8 10 
8 11 
8 12 
8 13 
8 Ui 
8 15 
8 16 
8 17 
8 18 
8 19 
8 20 
8 2! 
8 22 
8 23 
A 24 



1 0,0 

1 0,0 

1 0,0 

I 0,0 

1 0.0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0.0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 

1 0,0 



59,00 
57.00 
56.00 
56,00 
55,00 
55,00 
58,00 
61,00 
66.00 
68,00 
71.00 
7*. 00 
7a, 00 
76,0 
78.00 
77,00 
77,00 
78,00 
77,00 
75,00 
72,00 
68,00 
66,00 
64,0 



0,00040 
0.000*9 
0,00039 
0,00039 
0,00039 
0,00039 
0,00039 
0.00038 
0,00038 
0,00038 
0,00038 
0,0003* 
0,0003» 
0,00038 
0,00038 
0,00037 
0,00037 
0,00037 
0,00037 
0,00037 
0,00037 
0,00037 
0,00037 
0,00037 



4,86115 INCME3 



62.00 
61,00 
60,00 
58 00 
58,00 
56.00 
61,00 
64,00 
68,00 
71,00 
73.00 
74,00 
76 00 
77,00 
78,00 
79,00 
78,00 
78,00 
78,00 
74, 00 
71,00 
b9,00 
6 7,03 
65,0 



0,00037 
0,00037 
0,00036 
0,00036 
0,00036 
0,00036 
0,00036 
0,00036 
0,00036 
0,00036 
0,00036 
0,00036 
0,00035 
0,00035 
0,00035 
0,00035 
0,00035 
0,00035 
0,00035 
0,00035 
0,00035 
0,00035 
0,00034 
0,00034 



0,59? 
0.578 
0,56b 
0,554 
0.544 
0,533 
0,524 
0,515 
0.506 
0.498 
0,49 
0.483 
0.476 
0.469 
0,462 
0,456 
0,450 
0,44 4 
0,438 
0,432 
0.4 26 
0,420 
0,414 
0,408 



0,403 
0,397 
0,391 
0,385 
0.379 
0,373 
0,366 
0.360 
0.354 
0,347 
0.341 
0,334 
0,328 
0,321 
0.314 
0.307 
0,301 
0.294 
0.287 
0,280 
3.273 
0,2b6 
0.259 
0,253 




► , 20225 
1,19997 
0,19778 
0,19567 
0,19364 
0,19166 
0,18974 
0,1878b 
0,18601 
0,18419 
0,18239 
0,18059 
0,17881 
0,17702 
0,17523 
0,17344 
0,17163 



0,16981 
0,16798 
0,16612 
0,16425 
0,16235 
0.16043 
0,15850 
0,15654 
0,15457 
0,15258 
0,15059 
0,14858 
0,14658 
0,14457 
0,14256 
0, 1405b 
0.138S7 
0,13h58 
0.1346? 
0,13266 
0,13073 
0,12881 
0.12692 
0.12505 



8,11 
8,11 
8,12 
8,13 
8,13 
8,14 
8,15 
8,15 
8,16 
6,16 
8,17 
8,17 
8.1H 
8,18 
8, 19 
8,19 
8,20 
8,20 
8,21 
8,21 
8,22 
8,22 
8.23 
8,23 



8.24 
0,24 
0,25 
0,25 

0,26 
8,27 
8.27 
8,28 
8,28 
8,29 
8,30 
6,30 
8,31 
8,32 
8,33 
8,33 
8,34 
8.35 
8.3b 
8,36 
8,37 
8,38 
8,39 
8,40 



2,53 
2.55 
2,56 
2.58 
2,59 
2,60 
2,61 
2,62 
2,64 
2,65 
2,66 
2,67 
?,6H 
2,69 
2.70 
2,70 
2.71 
2.7? 
2,73 
2.74 
2.75 
2,76 
2,77 
2,78 



0,23614 
0.23182 
0,22779 
0,22403 
0,22052 
0,21725 
0,21417 
0,21129 
0,20857 
0.20600 
0,20355 
0,20123 
0,19900 
0,l9b87 
0,19481 
0,19281 
0,19086 
0,18897 
0,18711 
0,10527 
0,10346 
0,18167 
0,17988 
0,17810 



2,79 
2,80 
2.02 
2,83 

2,04 
2,05 
2,86 
2,88 
2,09 
2,90 
2,91 
2,93 
2,94 
2,96 
2,97 
2,99 
3,00 
3,02 
3.03 
3,05 
3,06 
3,08 
3.10 
3.11 



0,17632 
0,17454 
0,17274 
0,17094 
0,16913 
0,16730 
0,16545 
0,16359 
0,16171 
0,15982 
0,15790 
0,15598 
0,15405 
0,15210 
0,15015 
0,14820 
0,14b?4 
0,14429 
0,14a35 
0,14041 
0,13848 
0,13656 
0,13466 
0,13278 



0,0 
0.0 
0.0 
0,0 
0,0 
0.0 
0,0 
0,0 
0,0 
0.0 
0.0 
0.0 
0,0 
0,0 
0,0 
0,0 
0.0 
0,0 
0.0 
0,0 
0.0 
0.0 
0,0 
0,0 



0,0 

0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0.0 
0.0 
0.0 
0.0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 



1,79842 
II a 79680 
31,79910 
$1,79956 
■1,79994 

1,80032 

■ 1,80070 

■ 1,80109 
1,80147 
1,8016b 

1,80222 
1,80258 
1.H029S 
1.H0331 

■ 1.60360 
■1.80405 

■ 1,80441 
1,80478 

■ 1,80515 
1,80551 

■ 1,80588 

■ 1,80624 

■ 1,60661 
1,00698 



31,80734 
31,80771 
31 ,80806 
31,80841 
31,80876 
31,80911 
31,60946 
31,80981 
31,81017 
31,81052 
31,81087 
31,81122 
31,81157 
31,81192 
31,81227 
31,81262 
31,8129/ 
31,81331 
31,81364 
31,81396 
31,81432 
31,81465 
31 ,B149<J 
31.81532 



0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 



0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0.0 
0,0 
0,0 
0,0 
0.0 
0,0 
0,0 
0,0 
0,0 
0,0 



0,100 
0,100 
0.101 
0.101 
0.101 
0.101 
0.101 
0,101 
0,101 
0,101 
0.101 
0.101 
0.101 
0,101 
0,101 
0,101 
0,101 
0,101 
0,101 
0.101 
0,101 
0,101 
0,101 
0,101 



9fc 01 

9fc 01 

OE 01 

Ot 01 

OE 01 

IE 01 

IE 01 

2E 01 

2E 01 

?-t 01 

3£ 01 

3E 01 

HI ul 

4fc 1 

5fc 01 

5fc 01 

5fc 1 

6k 01 

6E Ol 

6E 01 

7k 1 

7k 01 

Of. 01 

0E 01 



0,1018k 
0.1019E 
0,10l9fc 
0,l0l9fc 
0.1020E 
Q.1020E 
0.1020E 
0.1021E 
0.1021E 
0,1022k 
0,102?E 
0. 1022k 
0,1022fc 
0.1023E 
0, 1023k 
0, 1023k 
0,1024fc 
0,1024E 
0,1024k 
0, 1025k 
0.1025E 
0,1025k 
0,1025k 
0,l02bfc 



0,68b6k 
0,6881k 
O,h09bk 
0.b908k 
0.69?lk 
0,6934k 
0,6946E 
0.6961F 
0,6974k 
0,b987fc 
0.7000E 
0.7013F 
0,702bk 
0,7039k 
0.7052F. 
0.,7ub«4k 
0,7077k 
0,7090k' 
0,71 Oik 
. 7 1 1 bk 
0.712BE 
0.7140E 
0.7153E 
0.7165k 



02 
Ot! 
02 



02 
02 
02 
02 
02 
02 
Od 
02 o 
02 
02 
02 

Od U 

02 
02 
02 
02 
02 
02 u 
02 
02 
02 



,8744k 00 

,«753fc 00 

,67b2E 00 

,87/Ot 00 

,0779* 00 

,8788k. 00 

,8797£ 00 

,6806E 00 

.8814E 00 

,68236 00 

,8832fc 00 

,8840k 00 

,6849t 00 

,«857k 00 

,B8b6fc 00 

,08f4fc oO 

,88H3£ 00 

,689lk 00 

,fl900E 00 

,89o8k 00 

,69lbE 00 

,69«>5k. 00 

.8933E 00 

,8941k 00 



01 0, 

01 0, 

01 0, 

01 0, 

01 0, 

01 0, 

01 0, 

01 o, 

01 0. 

01 0, 

01 0. 

01 o, 

Ol 0. 

01 o, 

01 o, 

01 o, 

01 0, 

01 0, 

01 0. 

01 0, 

ol 0, 

01 0, 

01 0, 

01 P, 



7170E 
7190F 
7202k 
7214E 
7227k 
7239k 
7251k 
7263k 
7275k 
72«7k 
7299k 
7311k 
7323k 
7334E 
754bk 
7358k 
7369k 
7381k 
7393k 
7«04k 
7«16E 

7438k 
745uE 



02 0, 
02 0, 
02 0, 
02 0, 
02 0, 
02 0. 
02 0, 
02 0, 
02 0. 
Orf w, 
02 0, 
02 0, 
02 0, 
02 0, 
02 0, 
02 0, 
02 0, 
02 0, 
02 0, 
02 0, 
02 u, 
02 0, 
02 0. 
02 0, 



89b0fc 00 

8958E 00 

0966E 

0974fc. 

0982b 

6990E 

6998E 00 

9006E 00 

9014fc. 

9022k 

9030fc 

9038t 

90«6k 

9054E 

9062k 

90b9fc 

90/7E 

9065E 

9092k 

910O£ 

9 1o8k. 

9115k 

91<i3t 

9130k 



00 

oo 
00 

00 



00 
00 
00 
00 

uo 

00 
00 
00 

00 
00 
on 
00 

00 

00 
00 
00 



OLOOKd, 8) a 



3.71568 INCHES 



APPENDIX C 
TVA MODEL 
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APPENDIX C - TVA Model 

MODEL PARAMETERS AND CONSTANTS 

Primary Model Parameters 

B = a volumetric parameter used to preserve mass balance 

A = a surface runoff volume parameter for winter storms 

D = a surface runoff volume parameter for summer storms 

GWK = a ground water volumetric parameter 

TDSRO= a routing parameter 

Secondary Model Parameters That May Be Optimized 

SROK = a surface runoff recession constant 

GWMIN= a threshold constant used in allocating ground water 

BHORP= soil B horizon permeability 

GROKW= an interflow or winter recession constant 

GROKS= a summer ground water routing contant 

Model Constants 

ACREIN= drainage area in square miles 
WCEPT = winter interception capacity 
SCEPT = summer interception capacity 
AHORD = moisture capacity in the soil A horizon 
GWDOR = a ground water reservoir allocation constant 
FALL, WINTER, SPRING, SUMMER = day of the year for the beginning of the 
respective season 
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